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We investigate the possible decay of protons in geodesic circular motion around neutral compact 
objects. Weak and strong decay rates and the associated emitted powers are calculated using a 
semi-classical approach. Our results are discussed with respect to distinct ones in the literature, 
which consider the decay of accelerated protons in electromagnetic fields. A number of consistency 
checks are presented along the paper. 



PACS numbers: 14.20.Dh, 95.30.Cq 



I. 



INTRODUCTION 



It is well known that according to the particle standard 
model inertial protons are stable. However, this is not so 
if the proton is under the influence of some external force 
because in this case the accelerating agent can provide 
the required extra energy, which allows the proton to 
decay. To the best of our present knowledge, the first 
ones to consider the decay of accelerated protons and 
similar processes as 



P 



P 7T 



(1.1) 



were Ginzburg and Zharkov I n Ref- Q the proton 

is described by a classical current with a well defined tra- 
jectory while the pion is field quantized. This approach 
is accurate in the no-recoil regime, i.e. when the relevant 
parameter \ = o,/m p involving the p +, s proper acceler- 
ation a and mass m p is less than unity. At the same 
time, Zharkov (see also Ref. for a recent review) 
investigated the process 



to investigate quantum processes in such a background. 
The study of particle processes in the presence of strong 
electromagnetic fields should be important to the analy- 
sis of certain aspects of high energy cosmic ray physics 
produced in pulsars and magnetars. In such intense mag- 
netic fields (H ~ 10 12 — 10 17 G) the strong coupling of 
protons and neutrons with mesons cangenerate p's and 
7r's with a non negligible intensity @-0- 

In the proper regime (i.e. where backreaction effects 
are not important), the reaction rate associated with pro- 
cesses (|l.lf) when the p + is in circular motion and i|1.2|) 
when it is under the influence of a magnetic field coincide. 
Notwithstanding, this is not so for the processes (|1.3(l - 
and 



p — ► n 7t , 

+ a + 

p —> n e v . 



(1.5) 
(1.6) 



respectively (where the p" 1 



are described by a classi- 



cal current). This is a consequence of the fact that in the 
classical current approach both, proton and neutron, are 



P 7T 



(1.2) 



(0 



and the strong and weak proton decays 



-{- M + 

p — > 77,71" , 

+ A * + 

p — > n e v , 



(1.3) 
(1.4) 



respectively, in the presence of an electromagnetic field 
Ap, where all particles are field quantized. For this pur- 
pose, it was used the comprehensive formalism developed 
by Nikishov and Ritus (see also || ) , which allows one 
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FIG. 1: The process p + — > p + % is represented in the pres- 
ence of background gravitational (h^) and electromagnetic 
(A 1 *) fields. In the proper regime, this process can be well 
described in both backgrounds by the classical current ap- 
proach, where the proton is assumed to have a well defined 
worldline with acceleration a and the pion is field quantized. 
Notice the similarity of the proton behavior in (i) and (ii). 
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FIG. 2: The process p + — > n n is represented in the pres- 
ence of background gravitational (h^) and electromagnetic 
(A 11 ) fields. Only in the first case it can be well described, in 
the proper regime, by the classical current approach labeled 
by the nucleons' proper acceleration a. Notice the difference 
of the proton-neutron behavior in (i) and (ii). 



usually assumed to follow the same trajectory in contrast 
with what really happens in the presence of a background 
magnetic field (notice the difference between Fig. ^ an d 
Figs. Hand El). 

This raises the question about what is the physical situ- 
ation simulated by the classical current method when ap- 
plied to the processes (|1.5H and Ijl.tifl . where one considers 
that only mesons and leptons are field quantized. Once 
protons and neutrons are described by a common current, 
one should look for a situation where they are mainly 
undistinguishablc. This is what happens in gravitational 
fields according to the equivalence principle. (For early 
and recent investigations on geodesic synchrotron radi- 
ation from classical currents see Refs. || and @, re- 
spectively.) As a consequence, processes (|1.5f) and (|1.6f) 
should represent fairly well the strong and weak conver- 
sion of protons into neutrons when they orbit chargeless 
compact objects provided that the back reaction on the 
neutron is negligible. This is what we are going to inves- 
tigate in this paper. 

In our procedure, we take into account the proton- 
neutron mass difference, by introducing a semiclassi- 
cal rather than classical current. We will be following 
Ref. [lfj, where a semiclassical current was successfully 
used to model the decay of linearly accelerated protons 
in the study of the Fulling-Davies-Unruh effect [TlJ . The 
current is "classical" in the sense that the proton-neutron 
is associated with a well defined worldline and "quan- 
tum" in the sense that it behaves as a two-level quantum 
system. (A simplified related calculation, where all parti- 
cles are treated as scalars can be found in Ref. U3-) The 
calculation is performed in Minkowski spacctime and the 
gravitational field is described by a Newtonian-like cen- 
tral force. 

There is also an important difference concerning the 
proton decay when it is under the influence of a grav- 
itational field rather than of a magnetic one, which is 
worthwhile to call the attention. The physical scale for 
the proton decay is given by its proper acceleration a. 
Because of the pion mass, process (|l.bfl dominates over 




FIG. 3: The process p + — > n° e + v is represented in the pres- 
ence of background gravitational (h^) and electromagnetic 
(A^) fields. Analogously to the previous figure, only the first 
case can be well described, in the proper regime, by the clas- 
sical current approach. 



process ljl.5|l in the region m e + A/z < a < m-n + A/z and 
because of the magnitude of the strong coupling constant, 
process (|1.5|l dominates over process (|1.6ll in the region 
a > to,,. + A/z, where A/z = m n — m p . Now, in the pres- 
ence of a magnetic field H , the proper acceleration of a 
proton in circular motion can be written as a = jeH/mp, 
where 7 = E/m p is the usual relativistic factor given by 
the ratio of the proton's energy and mass. In the region 
a > m n + A/z, the strong process dominates over electro- 
magnetic processes and energy degradation through pho- 
ton emission does not play any relevant role. This is not 
so, however, in the region m e +A/i < a < m^+Afi, where 
electromagnetic processes dominate over the weak one 
and much of the proton's energy E can be carried away 
by the photons, driving its acceleration below the thresh- 
old m e + A/i. (Recall that 7 is proportional to E.) The 
situation is quite different in a gravitational field. Assum- 
ing Minkowski space, the proper acceleration a = R£l 2 j 2 
of a proton in circular orbit with radius R and angular 
velocity Q around a compact object with mass M can be 
written as a = (GAfft 4 ) 1 / 3 /^ - (GAffJ) 2 / 3 ), where we 
have used the Newtonian gravity relation R 3 £l 2 = GM. 
Then, as the orbiting proton emits photons descending to 
a more internal orbit with larger f2, its proper accelera- 
tion tends to increase rather than to decrease, in contrast 
to the electromagnetic case. Whether or not a proton de- 
cays along its inspiraling trajectory will depend on the 
mass of the central object and other details, which will 
be discussed further. 

Now one may wonder how accurate can be our re- 
sults when applied to quite strong gravitational fields. 
In principle, an exact calculation would require that we 
take into account the spacetime curvature in the par- 
ticle field quantization. However, as it was shown in 
Refs. 13 and [14( the results obtained assuming a full 
curved spacetime ruled by Einstein equations and a flat 
background endowed with a Newtonian attraction force 
should not differ by more than 20%-30% up to the inner 
stable circular orbit of a static black hole. This is going 
to suffice for our purposes. 
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The paper is organized as follows: In Sec. [H] we present 
the semiclassical current formalism. In Sec. lIIII we evalu- 
ate the scalar and fermion emission rates for a uniformly 
swirling current. In Sec. II VI we evaluate the correspond- 
ing radiated powers. In Sec. 0we use the previous re- 
sults to analyze the decay of protons orbiting chargeless 
compact objects. Consistency checks for our formulas 
are presented. We dedicate Sec. El to our final discus- 
sions. We assume Minkowski spacetime with metric com- 
ponents rj^ — diag(+l, — 1, — 1, — 1) associated with the 
usual inertial coordinates (t, x) and adopt natural units 
c = % = 1 throughout this paper unless stated otherwise. 



II. SEMICLASSICAL CURRENT FORMALISM 



Let us consider the following class of processes 



Pi 
Pi 



P2 g , 
P2 fi h , 



(2.1) 
(2.2) 



where a scalar g or a fermion-antifermion pair /1-/2 are 
emitted as the particle pi evolves into pi. The g, fi, 
fi, pi and pi's rest masses are m, mi, mi, Mi and Mi, 
respectively. We will be interested here in cases where 
m,TOi,m2 <C Mi, Mi. The particle emission will be as- 
sumed not to change significantly the four- velocity of pi 
with respect to p\. This is called "no-recoil condition", 
which is verified when the momentum of the emitted par- 
ticles with respect to the instantaneously inertial rest 
frame lying at pi satisfies |k r f| <C Mi, M 2 - Because 
m, mi, mi <C Mi, Mi, this implies that the energy of each 
emitted particle satisfies w r f <§; Mi, M 2 . As the typical 
energy w r f of the emitted particl es is of the order of pi 's 
proper acceleration a = y f \a ll a tl \, the no-recoil condition 
can be recast in the frame independent form ^| 



a < Mi, Mi 



(2.3) 



The particles pi and pi will be seen as distinct en- 
ergy eigenstates \pi) and \pi), respectively, of a two-level 
system. The associated proper Hamiltonian Hq of the 
particle system satisfies, thus, 



H \pj)=M J \ Pj ) , j = 1,2 



(2.4) 



We shall describe our pointlike particle system pi-pi in 
the process (|2.1|l by the semiclassical scalar source 

j(x) = [q(T)/u°(T)]S 3 [x-x(T)] (2.5) 

and in the process 11'. 21 by the vector current 

f(x) = [q(r) u»(t)/u°(t)} <5 3 [x - x(r)] . (2.6) 

Here x^(t) is the classical world line parametrized by 
the proper time r associated with pi-pi, u m (t) = 
dx^/dr is the corresponding four- velocity, and q(r) = 

e lH ° T qoe~ lH ° T , where qo is a self-adjoint operator evolved 



by the one-parameter group of unitary operators U(r) = 

e -iH T 

The emitted scalar g in the process l|2.1(l is associated 
with a complex Klein-Gordon field 



(2.7) 



*(*)= / d d k a k ^(i) + ^- k '(!t) 



while the emitted fermions fi — fi in the process <|2.2|) 
are associated with the fermionic one 

a=± J 

(2.8) 

where i = 1,2 labels the two fermions. Here &k (a^)' 
bki<7 (&k i( r)' (c k ) and d^ iCr (d^ a ) are annihilation 
(creation) operators of scalars, fermions, antiscalars and 
antifermions, respectively, with three-momentum k = 
(k x , k y ,k z ) and energy lj = Vk 2 + m? for the scalar, and 
kj = {kf,kV,kf) and = ^/k 2 + mf for the fermions. 

< ^k ± "' ) ano - Wc^*' are positive [(+w), (+Wj)] and negative 
[(— ui), (— L>i)] frequency solutions of the Klein-Gordon 
(□ - m 2 )(f> ( ^ u) = and Dirac (i^d^ - mi)^^ = 
equations, respectively, where a labels the fermion polar- 
ization. 

Next, we minimally couple the fields to our semiclas- 
sical source (I2.5[) and current 1)2.6(1 according to the ac- 
tions d-Gili 



<?W = j d A x j(x) $(x) + & {x) 
for the scalar and 

S ( j f) = J d A x~3^ir{c v -c Al 5 )i>i 
+ #27 M (cv - cat 5 )*!) 



(2.9) 



(2.10) 



for the fermionic cases, (I2.1|l and (|2.2(l . respectively, 

where $ = and cy = ca = 1 in the processes 

here analyzed. 

The transition amplitudes at the tree level for the pro- 
cesses ()2.1|l and ()2.2ll are given by 



-4k = (Pi\®(g k \ S { j s) |0)<8|jji) , (2.11) 



and 



K\Z = (Pal ® (/ik lCTl! / 2k2(72 | SP |0) ® \ P i) , (2.12) 
respectively. The differential transition probabilities are 



d 3 k 



= IA| 2 



G 



( S )2 
cff 



dr / dr' exp{iA/i(r — t) 

-00 J — 00 

(2.13) 



2{2nfu J_ 
+ik x [x{r) - x(t')]a} 
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and 
dVf^ P2 

1 z a 1 =±<7 2 =± 



o-2 |2 
k 2 l 



= ro ^ / dT / dT ' exp{*A A1 (r - r') 
+ i(fc 1 + fc 2 ) A [ a; (T)-x(r')]A} 

X {2k^k'^U^T)u L ,(T / ) - k°k 2a U»(T) Uf ,(T') 

+ te^ af3 k la k 2 pu^T)u„(T')} , (2.14) 

accordingly, where g^ Ql/ ^ is the totally skew-symmetric 
Levi-Civita pseudo-tensor (with 6 0123 = —1), k^k v 2 = 
(fcffc^ + Kk 2 l )/2, Afi = M 2 - Mi and G { c f if) = 
I (P2\qo\pi) | are the effective coupling constants for the 
scalar (s) and fcrmionic (/) channels. 



III. 



EMISSION RATES 



The world line of a particle with uniform circular mo- 
tion with radius R and angular velocity f2 as defined by 
laboratory observers at rest in an inertial frame with co- 
ordinates (t, x), is 



x»(t) = (t , Rcos(nt) , Rsm(nt) , 0) 



(3.1) 



and the corresponding four-velocity is 

u»(t) =7(1 , -RClsia(Qt) , RClcos(Qt) , 0) , (3.2) 

where 7 = (1 — i? 2 17 2 ) -1 / 2 = const is the Lorentz factor 
(v = RD, < 1), t = 7T, and a = x /-a (l a' 1 = Rfl 2 ^ 2 is 
the proper acceleration. Let us calculate now separately 
the scalar and fermionic emission rates associated with 
processes l|2.1|l and ll'2.2t . respectively. 

A. Scalar case 

First, let us analyze the process (|2.1|) . In order to 
decouple the integrals in Eq. (|2.13|l . we define new coor- 
dinates, 

cr=7(r-r')/2 and s = 7 (r + r')/2 , (3.3) 
and perform the change in the momentum variable 

fc"h-> fc"= (£,k) , (3.4) 

where 

UJ = UJ , 

k x = k x cos(fts) + k v sin(fts) , 
k y = -k x sin(f2s) + k v cos(fis) , 
\\j — - k ^ , 



which consists of a rotation by an angle Qs around the 
k z axis. Hence, we obtain from Eq. (|2.13l) the follow- 
ing transition rate per momentum-space element for the 
emitted scalar: 



d 3 k 



G 



(s)2 
off 



(2tt) 3 27 2 w 



/+oo 
da exp[i(A/xcr/7 
-OO 



(3.5) 



where TZ^ 1 ^ 2 = d'P Pl ~* P2 /ds is the transition probabil- 
ity per laboratory time and 



J"((7) = (a, , 2i?sin(ftcr/2) , 0) 
In order to calculate the transition rate 



d d k 



dKP 1 ^ 2 
d 3 k 



(3.6) 
(3.7) 



we use Eq. (|3.5|l and obtain 



G 



(s)2 
off 



2 7 2 (2tt) 3 J_ 



4-oo 



da e l Af * Gh I{a) , (3. 



where 



1(a) = / d 3 k 



j k x x x 



(3.9) 



and uj = \/k 2 + to 2 . In order to integrate Eq. (|3.9|l . 
we introduce spherical coordinates in the momenta space 
(k e R+,0 G [O,7r],0 S [0,2tt)), where fc* = fcsinflcos^, 
k y = k sin 9 sin 0, and k z — kcosO. By doing so, we 
obtain 

4tt '•+ 00 



due tuA sin 



Vw 2 - m 2 |X| 



where |X| = V^ATX* . Next, by redefining the fre- 
quency variable as uj = to cosh £, we obtain 

7( CT ) = z^p y +00 ^ e ™(x°cosh« + ixisi„h« )sinhe 

Now, we perform the change of variable ^ i— * 77 = e^, 
leading to 



/(a) 



drj(rj 



1) exp 



(y° + |Y|)77 



im(y° - |Y|) 



2r? 



(3.10) 



where we have introduced a small positive regulator e > 
in the integral as follows: 



X» h-> = (X° + ie, X\X 2 , X 3 ) 



(3.11) 



(Note that |Re(F°)| = > |X| = |Y|.) Then^v 
using expressions (3.471.11) and (8.484.1) of Ref. [l]f, 
we obtain 

, . — 2 7r 2 i msign(cr) cii / , , / \ 

/(a) = * y J gf (sign(a)TO V ^F^) , 

(3.12) 



where H[ ' (z) is the Hankel function of the first kind. As 
a result, by making the variable change a i— ► A = — acr/7 
and by defining Z 11 = (afj)Y 11 , the transition rate (|3.8|) 
can be cast in the form 



-iG { e ffh 2 a r-*o 



8nj 



dX e~ l Afl 



^\z) 



_ (3.13) 
where we have defined m = m/a, A/i = A/i/a, e' = 
ae/7 -C 1, z = — m7sign(A)\/ Z\Z X , and where 

Z» = (-A + ie', 0, -(2Ra/j) sinfOVy/^a), 0). (3.14) 

Eq. I|3.13|l is our general formula for the transition rate 
per laboratory time. 

In the physically interesting regime, where fh <C 1 it 
can be integrated using the following expansion for the 
Hankel function [l8|: 

H^\z) ps -— + 0(zlnz) for \z\ « 1 . (3.15) 

We note that for large enough |A|, \z\ > 1, Eq. (|3.15(l 
ceases to be a good approximation. [For instance, for 
7 2 > l/m > 1, we have that |z| > 1 for |A| > 1/\/12to, 
while for l/m ^> 7 2 3> 1, we have that \z\ > 1 for 
|A| > 1/(7™).] Notwithstanding, this is not important 
because the error committed in this region is small to 
affect the final result provided that iii C 1. Hence we 
write Eq. Ij3.13|l for fh -C 1 in the form 



TZ Pl 



-^cff a I dX e 



47T 2 7 3 



{Z X Z X ) 



(3.16) 



where 



Z X Z X = (A — it') 2 - {2Rah) 2 sin 2 (f7A7/2a) . (3.17) 

In order to solve this integral, we expand Z X Z\ for rel- 
ativistic swirling particles [13, i.e., 7 » 1 (recall that 
R = v 2 j 2 /a, = a/(vj 2 ), and v = \/T— 7 -2 ): 

Z X Z X ps T ^(A + ^\/3A + )(A + ^\/3yl_)(A-^\/3B + ) 

(3.18) 



where 



:(A-n/3B_) , 



AT = lTV 1 + 2? /^ / 3 



and 



BTEl T yi-2?/\/3 

with 1. For |A| > 2vj, where the expansion ceases to 
be a good approximation, the integral contributes very 
little again and, thus, will not have any major influence 



c 



FIG. 4: The integration path in complex plane. The Xa±,b± 
refers to the poles for each of the terms in parentheses of the 
Eq. iSTTHI . 



in the final result. Thus, the integral in Eq. I|3.16|l can 
be rewritten in the complex plane: 



TZ Pl 



r<(s)2 

eff a I dX e 



47T 2 7 



2~,3 



-i A/i A 



z x z x 



(3.19) 



where the complex integration path C is clockwise ori- 
ented (C = (-L,L) U {Li f e 7 9 e [-tt,0]},L -> 00) as 
shown in Fig. 0] Eq. H3.19(l can be performed, then, by 
using Cauchy's residue theorem leading to 



K Pl 



^eff a -2V3A~I* 



8V37T7 

which is valid for fh <C 1 and 7> 1. 



(3.20) 



B. Fermionic case 



Next, let us compute the transition rate associated 
with the fermionic fa — fa emission of process (|2.2|l . Af- 
ter performing the variable changes 1)3. 3|) and 1)3. 4JI . one 
obtains from Eq. (|2.14|) . 



dK Pl ^ P2 



2G. 



(/)2 
eff 



/+00 
<i<7 exp[z(A/icr/7 
-00 



d 3 kid 3 k 2 (27r)%iS 2 ^^oo 

+ k 2 ) p X fl (a))][R 2 il 2 (^ 2 ~k(k z 2 )cos{il(j) 
-R 2 n 2 (kfk$ - fcffcf) + (wiC 2 + ki ■ k a ) 

- 2 #fi (Sifcf + u 2 kl) cos( Qo-/2 ) 
+ 2iRfl(k 1 x k 2 ) :E sin(f7 ( T/2) 

- i i? 2 ft 2 - £ 2 %) sin( fla )] , (3.21) 

which is the laboratory transition rate per momentum- 
space element associated with each emitted fermion and 
X p is given in Eq. (|3.6|l . By integrating over all mo- 
menta, the transition rate can be rewritten in a more 
convenient form as 



2 G 



Cf)2 r +oo 

eff 



(27T) 



da e l A " Gh G^A pu (3.22) 



6 



with 



G 



{lis 



dh dh 
dX^ dX v 



(3.23) 



where the index I — 1, 2 is used to distinguish the fermion 
in the final state to which we are referring and lo\ — 

k 2 + to 2 . Also 



and 



d 3 k, 



(3.24) 



1 + R 2 fl 2 cos(ftcr) 

o i - R 2 n 2 

-2Rflcos(fla/2) 
iR 2 fl 2 sin(f)cr) 



-2M}cos(ficr/2) -iR 2 VL 2 sin(ficr) 



1 + i? 2 ^ 2 2iRfisin(Ocr/2) 

-2iRtt sin(ficr/2) 1 - R 2 tt 2 cos(ftcr) 



(3.25) 



In order to compute Eq. (|3.24|) . we introduce spherical coordinates in the momenta space (fc; S R + ,0i e [0, 7r],0/ S 
[0, 27r)), where fcf = ki sin0; cos(/>;, fc^ = fc; sin6*; sin</>;, and fcf = ki cos 9i and perform the same steps of the previous 
section which led Eq. (|3.9|l into Eq. I|3.12[) . We obtain, thus, 



_ -2n 2 imi sign(g-) ^(i) 



HP (sign(a)m,Vn^ r ) 



where is dehncd in Eq. (|3 . 1 1|) . By introducing again Z^ = (a/~f)Y^ and the variable a 
transition rate (|3.22|l can be cast in the form (see also expression 8.472.4 in Ref. |l8j) 



(3.26) 
-aa/j, the 



_ |feff m l m 2 a 7 



dX e - iA ^Z»Z»A lw ^-^. 



z l 



(3.27) 



where to,/ = mi /a, Au = A/i/a, e' = ae/7 <C 1, 2; = — TO/7sign(A)\/ Z\Z X and is given in Eq. (|3.14|l with 



.4, 



" 1 + R 2 n 2 cos{tl~f\/a) -2i?ficos(Q7A/2a) iR 2 n 2 sin(0 7 A/a) 

l-i? 2 ft 2 

-2M7cos(ft 7 A/2a) 1 + R 2 Q 2 -2ii?ftsin(ft 7 A/2a) 

-iE 2 Q 2 sin(0 7 A/a) 2iRtt sin(ft 7 A/2a) 1 - i? 2 ^ 2 cos(0 7 A/a) 



(3.28) 



r 



This is our general expression for the laboratory reaction 
rate associated with the process (|2.2|) . 

Next, we cast Eq. ()3.27(l in a simpler form in the regime 
where to/ <C 1. For this purpose we use the expansion 
for the Hankel function 

H ( 2 ] {zi) « —4 - - + C(z; 2 lnz ; ) for |z/| < 1 (3.29) 

7TZ 7 7T 



the scalar case, the integral above is performed in the 
complex plane along the path given in Fig. ^ 



dXe -i-A,X Z " ZVA ^) 



8tt 4 7 7 c "" C (Z*Z X ) 
16 4(m 2 + m|) 

1 \Z X Z*) 2 + j 2 Z x Z* 



(3.31) 



and use a similar reasoning presented below Eq. (|3.15() 
[with the suitable identifications z — > 2;; and to, — > to,/] to 
obtain 



T cff 

Btt 4 7 



{z x z x y 



16 



4(to 2 + TO?,) 



(3.30) 



,7 4 (^a^) 2 7 2 ^ A 

where Z Z x is given in Eq. {nUBt for 7 > 1 (recall that 
i? = v 2 j 2 /a, Q = a/(v"f 2 ), and u = \/l — 7~ 2 ). As in 



Then, by using the Cauchy's residue theorem, we obtain 
for in 1, fri2 1 and 7 ^> 1 



G ( J s 'a 5 exp(-2V3A^) 
1728tt 3 7 



[49\/3 + 102A^ 



+30\/3A/i 2 + 12 A// - 39 a/3 (to 2 + m|) 
-90 A/i (to 2 + rh\) - 36 V3 A7t 2 (to 2 + m|) 



(3.32) 



This is easy to note that Eq. I|3.32|) is positive definite 
and decreases as rhi and TO2 increase. It is not difficult 
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to show, as well, that it also decreases as A/j, increases, 
as expected. 

As a check of our approach, let us use our formulas 
to analyze the usual /3-decay: n° — > p + e~ v. The mean 
proper lifetime of inertial neutrons is 887 s |2fJ|. Thus, 

n- 
f 



TC^ P = n n f ~" p {n -> 0) = ft/887 s leads to 



K': 



5.46 x 10" 



Gp MeV 5 



(3.33) 



where Gp = 1.166 x 10 -5 GeV~ 2 is the Fermi cou- 
pling constant poj . Clearly, we cannot use our expres- 
sion H3.30fl to calculate the reaction rate of the /3-decay, 
since it is not valid for inertial neutrons. However, 
can be derived in this case directly from Eq. H3.22J) by 
making £1 = in Eq. (|3.21|l . This is achieved by a change 
of the momentum variables as shown in Eq. (|3.4() . After 
performing the corresponding integrations in the angular 
coordinates and in ui ei we obtain 



Gy 

IT 



A/i— m e 



dio v u>l {Afj, - oj v ) 



x \J (An-aj u y 



(3.34) 



where we have assumed massless neutrinos, m v — 0, and 
G pn = Ggg ■ By evaluating numerically Eq. (|3.34() with 
m e = 0.511 MeV and Afi = (m n - m p ) = 1.29 MeV, we 
obtain 



K 



n — >p 



1.81 x 10" 



G 2 pn MeV 



This is to be compared with Eq. (|3.33|) . where G pn is 
to be identified with Gf- The reason why both results 
are not identical can be traced back to the fact that 
the nucleons are treated here semiclassically and have 
only approximately the same kinetic energy content: the 
no-recoil condition only models approximately the real 
physical situation. Notwithstanding, this suffices for our 
present purposes. 



IV. EMITTED POWER 
A. Scalar case 

Next, we calculate the radiated power 
W? 1 ^ 2 = / <i 3 ki ' 



d 3 k 



(4.1) 



associated with the emitted scalars as measured in the 
laboratory frame. Eq. I|4.1|l can be rewritten as 

= eff / d(J e z An ah J( a \ (4.2) 

(27T) J 7J-oo 



where 



= / d A ke 



3w k A X A 



(4.3) 



and is given in Eq. (|3.6|) . In order to integrate J (a), 
we follow closely the approach, which drove Eq. (|3.9|) into 
Eq. I-TT21 : 



J(a) 



2TT'm z Y (1) ( . 



(sign(a)mV^ r ) , (4.4) 



where V is given in Eq (|3.11|l . Now, by introducing 
again a h-> A = —aa/j and = (a/^Y^, Eq. (|4.2|l can 
be cast in the form 



87T 2 



(4.5) 

where z can be found below Eq. Ij3.13|l and is given 
in Eq. Q3.14p . Eq. (|4.5|) is the general expression for the 
radiated power associated with the emitted scalars. 

The expression above can be simplified in the limit 
m 1. For this purpose we use the expansion (see 
Ref. [3) 



4i 



(4.6) 



for \z\ < 1. Then, by letting Eq. (|4*1)|) in Eq. j|4~5|l . we 

can perform the remaining integral in the complex plane 
along the path of Fig. 0] to obtain the emitted power in 
the regime rh <C 1 and 7 ^> 1 

W?^ » ( 1 + ^AA . (4.7) 

This is in agreement with the expression obtained by 
Ginzburg and Zharkov Q (see also Ref. 0) in the due 
limit, i.e., A[i — ► 0. 



B. Fermionic case 

Further, we calculate the radiated power as measured 
by observers at rest in the laboratory frame associated 
with each fermion I — 1,2: 



fW J J d3 kld 3 k2 

Eq. I|4.8|l can be rewritten as 



(4.8) 



9 /~<(/) 2 r+cc 

— ^L. / dae'^^ff^, (4.9) 
( 27T ) J-00 



W f(l) ~ 



where we have chosen (with no loss of generality) 1 = 1, 
i.e., we are computing the radiated power associated with 
the fermion with mass mi. Here 



where 



dJi dl 2 
dX» dX v 



Ji(a) = I d 3 k ie lkl x > 



(4.10) 



(4.11) 
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and I2 is given in Eq. (|3.26() with 1 = 2. The result of 
Eq. (HE) : 



2ir 2 m 2 Y (i) 



(sign( ( r)m 1 0^ r ) (4.12) 



is obtained by inspection after comparing Eq. I|4.3|) with 
Eq. lj4~TTj) and Eq. (|4~4*|) with Eq. lj£T2)l . respectively, 
where 2;; is defined below Eq. 1(3. 27J) . and and are 
given by Eqs. I|3.14fl and (|3.28() . respectively. By letting 
Eqs. ljrT2l)jl (with I = 2) and gJ3 in Eq. (jCTIjl . we 
rewrite the emitted power (|4.9|) in the form 



W Pl^P2 

w fW 



8tt 2 7 - 2 

(1) 



dX e 

■00 



>1 



(4.13) 



This is our general formula for the total emitted power 
associated with the fermion 1 = 1. 

In the limit fhi <C 1, we can rewrite W P }^ P2 by using 



the expansions l|3.29|l and (see Ref. |l8j1 
2i zii 



H 3 ^)~-— f 



TTZl 



An 



+ 0{zf\nzi) 1 (4.14) 



for \ zi\ < 1. Thus, by letting Eqs. l|3~2l?j) and in 
Eq. (|4.13(l . we can perform the remaining integral in the 
complex plane along the same path shown in Fig. 01 and 
obtain the emitted power for fhi, 777.2 <C 1 and 7> 1: 



W Pl ^ P2 ■ 
W fW 



3456tt 3 



320 + 241\/3 A/j 



+246A/J + 46^3 A^ + 12A/Lt - 48(m 2 + 57772) 
-3\/3A/j(17m? + 657772) ~ lSA^ 2 ^ 2 + 13m 2 ,) 
-24V3A[i 3 (ml + 2ml)} ■ ( 4 -15) 



Clearly, W P ^ P2 is obtained by exchanging mi < — ► 7772 
in Eq. I|4.15|l . This is important to note that Eq. (|4.15(l 
is positive definite and decreases as tth, 777,2 and Afi in- 
crease, as expected. 

As a consistency check of our Eq. I|4.15|l , let us apply 
it to analyze the emission of neutrino- antineutrino pairs 
from accelerated electrons: — > e~ v e v e and compare 
the results in the proper limit with the ones in the lit- 
erature obtained when the electrons are qua ntized in a 
background magnetic field (see, e.g., [2l|-|25j and refer- 
ences therein). (Comprehensive accounts on 7 and v — v 
synchrotron radiation emitted from electrons in magnetic 
fields can be found, e.g., in Ref. |2al and Sec. 6.1 of 
Ref. [27j , respectively, and in Ref. |28j.1 The fact that we 
are assuming that our source is under the influence of a 



gravitational force rather than being immersed in an elec- 
tromagnetic field is not relevant in this particular case, 
since the neutrinos are chargeless (An account on the 
degradation of the neutrinos' energy in strong magnetic 
fields can be found in Ref. [29|.). As a consequence, our 
results and the aforementioned ones in the literature are 
expected to be in good agreement in the no-recoil regime 
X = a/m e <C 1. The total radiated power of neutrino- 
antineutrino pairs from circularly moving electrons in 
a constant magnetic field B with proper acceleration 
a = ^eB/m e -C m e (no-recoil condition) can be easily 
calculated W^ p = 5 (2 C v + 23 C\)G 2 F mix 6 / '(lOftrfj 
from the differential emission rate given, ejj;., in Ref. |25j 
or Ref. [13. Then (see Eq. (6.6) in Ref. \lj$), 



W^ p = 1.1 x 10" 



G 



where we have used that the vector and axial contribu- 
tions to the electric current are Cy = 0.93 and C\ = 
0.25 |2bT |. respectively, and x = a/m e <C 1. This is to 

be compared with the result obtained from Eq. (|4.15|) by 

( t\ 

defining G e „ = G cS and assuming A/i = m v = 0: 
^«lx 10- 2 G^ a 6 , 

where G e „ is the corresponding effective coupling con- 
stant, which is to be associated with the Fermi constant. 



V. PROTON DECAY 

Now, let us use our results to analyze the weak and 
strong proton decay processes l|1.6l) and (|1.5|) . respec- 
tively. Our formulas l|X77jl and ij4~T3|) . and (|XT3|) 
and l|4.5|) associated with the weak and strong reactions, 
respectively, are quite general although cumbersome to 
compute. Happily, we can use the much more friendly 
ones: <|3~32|) and (|4~T5j) . and (|3~2U|) and l|3"7)l. which 
are valid in the physical regime where processes <|1.6|) 
and l|1.5|) are more important. In the region where 



777 e «a« m,r 



(5.1) 



with TOtt being the tt + mass, the reaction (|1.6|1 has a non- 
negligible rate and dominates over the reaction i|1.5|) . In 
this case, Eqs. <|3.32|l and (|4.15|l can be used provided 
that 7>1. Now, in the region where 



m v <C a <C 777 p , 



(5.2) 



the reaction l|1.6|) is overcome by the strong one l|1.5fl . in 
which case Eqs. (13.20(1 and (|4.7|l should be used. Next, 
we look for orbits around compact object, where condi- 
tions l|5.1|) and l|5.2|l are verified. 

Let us begin by rewriting the proper acceleration of 
the proton in Minkowski space, a = R£l 2 j 2 [see below 
Eq. (|3.2[l ]. in the form 



r 



1 . 



(5.3) 
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log[a/(l MeV)] 

FIG. 5: The proton mean proper lifetime r associated with 
process II, til is plotted as a function of its proper acceleration 
a, where m e < a < m w . 




log[a/(l MeV)] 

FIG. 6: The proton mean proper lifetime r associated with 
process 11.51 is plotted as a function of its proper acceleration 
a, where m w < a < m p . 



Now, we use General Relativity to obtain the proton's 
energy per mass E/m p and angular velocity d4>/dr s as 
calculated by a static observer lying at rest at the same 
radius of the particle orbit around a compact object with 
mass M. E/m p and d<fi/dT s are to be identified with 7 
and fl in Eq. I|5.3[l . respectively, to obtain the proper ac- 
celeration a. Once we have a and 7, we use Eqs. <|3.32ll . 
(|4.15|l . (|3.2U|) and l|4.7|l to calculate the relevant decay 
rates and emitted powers. The results obtained in this 
way should be associated with the values defined by the 
static observers at the radius of the particle orbit. These 
ones differ from the reaction rates and emitted powers as 
measured at infinity by red-shift factors. In order to ob- 
tain (i) the reaction rates and (ii) the emitted powers at 
infinity from the ones measured by the static observers at 
the radius of the particle orbit, one should multiply the 
latter ones by (i) yjl - 2GM/r s and (ii) l-2GM/r s , re- 
spectively. Although we can only capture with this pro- 
cedure part of the influence of the spacetime curvature, 
its suitability as an approximate approach is justified by 
comparing the results which it provides with the ones ob- 
tained with full curved spacetime calculationSjWherever 
the latter ones are available, as, e.g., in Ref. |13| . 

The line element external to a spherically symmetric 
static object with mass M, which includes Schwarzschild 
black holes, can be written as 



and 



dS 2 = 1 - 



2GM\ 



dt 2 - 1 - 



2GM\ 1 



r I 



dr 2 - dY 2 , 



where g?£^ = r 2 (d9 2 + (sin 8) 2 dcj) 2 ) . According to Gen- 
eral Relativity |3(j , asymptotic observers associate an an- 
gular velocity d(j)/dt a . 0i = \J GM / rf and an energy per 
mass ratio E a , Jm = (1 - 2GM/r s )/y/l - 3GM/r s for 
particles in circular geodesies at r = r s . Thus, static ob- 
servers at r = r s (9, <f> = const) associate the following 
corresponding values: 



E/mp = y/1 - 2GM/r s /^l - 3GM/r s . 
By letting dcj>/dT s — * and E/m p — > 7, we obtain 



7 = 



' 1 - 2GM/r s 
1 - 3GM/r s 



and [see Eq. 1)5. 3fl ] 



GM 
r 2 (l-3GM/r s ) 



(5.4) 



(5.5) 



d<p/dr s = v/GM/rf/Vl - 2GM/r s 



which will be used to evaluate Eqs. I|3.32fl . H4.15JI . (|3.2U|) 
and H4.7fl . whenever 7 » 1. We note that Eqs. I|5.4|) 
and (|5.5|) are monotonic functions, which approximate 
the correct values asymptotically and diverge at r s = 
3GM. This is so because according to General Relativ- 
ity, circular geodesic orbits at r s 1=3 3GM approximate 
lightlike worldlines. 

At the last stable circular orbit, r s = 6GM, we obtain 
from Eq. I|5.5|l that 

a/m e = 3 x 10- 16 (M G /M) . 

Thus, protons around black holes in stable circular or- 
bits 6GM < r s < 00 are not likely to decay unless the 
compact object is a mini black hole with the mass of a 
mountain: M < 10 17 g [see Eq. (|5~T|) ]. The fact that 
the smaller the black hole the more likely that protons 
decay at a fixed r s /{GM) is related with the fact that 
the smaller the black hole the larger the spacetime cur- 
vature, i.e. "gravitational field", at the same r s /(GM). 

In order to explore more realistic cases, where black 
holes have some solar masses, we have to consider protons 
at inner circular orbits, 3GM < r s < 6GM, which are 
unstable. By defining r s = 3GA/(1 + S) with 6 <C 1 to 
monitor how far from the most internal circular orbit (at 
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r = 3GM) the proton is, we rewrite Eqs. (|5.4() and (|5.5() 

as 

7«1/a/3£ (5.6) 

and 

a m 1/(9GM<5) . (5.7) 

By using Eqs. I|5.6[) and (|5.7|) in Eqs. (|5.1|l and (|5.2|) we 

obtain 

3 x 10" 19 (Af Q /M) < 5 < 9 x 10" 17 (M Q /M) (5.8) 

and 

5 x 1Q- 20 (M G /M) < 5 < 3 x 1O~ 19 (M /M) , (5.9) 

which are the intervals where the weak and strong pro- 
cesses would be favored, respectively. Thus, free protons 
in circular orbits around stellar mass black holes are likely 
to decay only if they are extremely close to the most in- 
ternal circular geodesic and stay there for long enough to 
decay. 

In Fig. 03 we plot from Eq. (|3.32J1 the proton mean 
proper lifetime r(a) — 1/T p ^ n associated with the pro- 
cess (|1.6|) . where T p f~ yn = jTZ^ n is the weak tran- 
sition probability per proper time and we have iden- 
( f ) 

titled G cS — G pn with the Fermi coupling constant 
G F = 1.166 x 10" 5 GeV -2 . We have plotted the proper 
lifetime t(o) rather than the laboratory lifetime t(a) in 
order to make it easier the comparison of this figure with 
Fig. 1 in Ref. 15]. In Fig. H we plot from Eq. (|3~27)|) 
the proton mean proper lifetime r(a) = 1/TJ - *™ associ- 
ated with process (|1.5fl . where rf^ n = ^W ! ~' n is the 

(s) 

strong transition probability per proper time. Here G c ^ 
is identified with the pion-nucleon-nucleon strong cou- 
pling constant g^NN, which is written in the Heavisidc- 
Lorentz system as \J g^ NN / (47r) ~ \/l4 (see, e.g., pj 
and Finally in Figs. \7\ and |H1 we plot the emitted 

power in the form of electrons and neutrinos as calculated 
from Eq. I|4.15[) and in the form of pions as calculated 
from H4.7[) . respectively. 

VI. DISCUSSION 

The decay of accelerated protons has attracted interest 
for long time. Astrophysics seems to provide suitable 
conditions for the observation of the decay of accelerated 
protons. Cosmic ray protons with energy E = "/m p 
3 x 10 14 eV under the influence of a magnetic field H ss 
10 14 Gauss of a pulsar have proper accelerations of an = 
jeH/rrip ss 200 MeV > m^. For these values of E and H, 
the proton are confined in a cylinder with typical radius 
R fa 7 2 /a_ff «2 x 10~ 2 cm <C Ih, where Ir is the typical 
size of the magnetic field region. Under such conditions, 
protons could rapidly decay through strong interaction 
before they lose most of their energy via electromagnetic 
synchrotron radiation 0. 
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0.5 1 1.5 2 

log[a/(l MeV)] 

FIG. 7: W e + and W v associated with process (11.611 are plotted 
as functions of the proton proper acceleration a with solid and 
dashed lines, respectively. 




2.2 2.4 2.6 2.8 3 

log[a/(l MeV)] 



FIG. 8: Wk associated with process 111.511 is plotted as a 
function of the proton proper acceleration a. 

Here we have considered the possible weak and strong 
proton decays under the influence of background gravi- 
tational fields. Reaction rates and emitted powers were 
calculated. We have concluded that they are unlikely 
to decay unless they orbit mini-black holes or they are 
pushed to highly relativistic geodesic circular orbits (and 
stay there for long enough to decay). This raises the ques- 
tion whether there would exist other astrophysical sites, 
where the decay rate could be larger. Perhaps the consid- 
eration of protons grazing the event horizon of black holes 
or entering properly the ergosphere of Kerr black holes 
extracting rotational energy from it would be worthwhile 
to be investigated. Notwithstanding because these cases 
would involve more complicated "trajectories" in a gen- 
uine general relativistic context, full quantum field the- 
ory in curved spacetime computations, rather than our 
semiclassical ones, would be desirable to provide more 
comprehensive results. 
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